Conservation and entanglement of Hermite-Gaussian modes in parametric 

down-conversion 
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We show that the transfer of the angular spectrum of the pump beam to the two-photon state 
in spontaneous parametric down-conversion enables the generation of entangled Hermite-Gaussian 
modes. We derive an analytical expression for the two-photon state in terms of these modes and show 
that there are restrictions on both the parity and order of the down-converted Hermite-Gaussian 
fields. Using these results, we show that the two-photon state is indeed entangled in Hermite- 
Gaussian modes. We propose experimental methods of creating maximally-entangled Bell states 
and non-maximally entangled pure states of first order Hermite-Gaussian modes. 

PACS numbers: 42.50.Dv, 03.67.Mn 



I. INTRODUCTION 

Recently, a great deal of attention has been paid 
to the higher-order Gaussian modes of the electromag- 
netic field. In the paraxial approximation, two interest- 
ing cases are the Laguerre-Gaussian (LG) and Hermite- 
Gaussian (HG) modes. These modes are solutions of the 
paraxial Helmholtz equation EI and are eigenstates of 
the free-space propagator. It has been shown that the 
LG modes carry orbital angular momentum in the form 
of an azimuthal phase e %i ^ in the transverse plane 0, 13 • 
Allen et al. 0] have shown that field modes with this 
type of phase dependence carry an orbital angular mo- 
mentum of £H per photon, where £ is the azimuthal beam 
index. 

These higher-order modes are of great interest in quan- 
tum information schemes, since they can be used to 
represent discrete 2?-state qudits. For example, the or- 
bital angular momentum of single photons in LG modes 
provides a possible qudit encoding scheme. The quan- 
tum number £ can be coherently raised or lowered us- 
ing holographic masks 0. One can measure the orbital 
angular momentum (LG modes) of single photons up 
to a desired accuracy using interferometric techniques 
[jj |(|- In, addition, Mair et al. have shown ex- 
perimentally that with spontaneous parametric down- 
conversion (SPDC) it is possible to create photon pairs 
entangled in orbital ang ular momentum and other theo- 
retical H H ES [H El [H El and experimental works 
E3 have followed, including the generation of entangled 
3-state qutrits El> Ey| • 

The HG modes may also be of use in quantum infor- 
mation schemes. In particular, the first-order HG and 
LG modes can be described and manipulated in a way 
that is analogous to linear and circular polarization of 
the electromagnetic field E3- Devices that act on the 



first-order transverse mode in a manner equivalent to 
polarizing beam splitters, half- wave plates and quarter- 
wave plates can be constructed using asymmetric inter- 
ferometers |l8l Il9j , Dove prisms |Tt| and mode convert- 
ers [23. The first-order modes can thus be used to de- 
fine qubits and the above devices implement single qubit 
rotations. Recently, Langford et al. have produced pho- 
tons entangled in first-order HG mode and performed 
quantum state tomography using holographic masks and 
single mode fibers [la]- Since the HG modes form an 
infinite-dimensional orthonormal basis, they too might 
be used to encode higher-dimensional qudits. 

Here we provide a theoretical description of the gen- 
eration of entangled HG modes for arbitrary HG pump 
beams. We show that there are restrictions on the parity 
and order of the down-converted HG fields. We intro- 
duce our notation in section II Al and briefly review the 
two-photon quantum state generated by SPDC in section 
ITT1 Our main results concerning the generation of corre- 
lated HG modes using SPDC are derived in section ITTT1 
including a general expression for the probability ampli- 
tude to generate combinations of different HG modes. In 
section IIVI we provide a proof that the down-converted 
HG modes are indeed entangled and we discuss the ex- 
perimental generation of Bell-states and non-maximally 
entangled pure states. 



A. HG modes 

For convenience, we will adopt the notation used in 
Ref. [2(J ■ The Hermite-Gaussian modes are given by the 
complex field amplitude 
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where the coefficients C nm are given by 

and H n {x) is the n th -order Hermite polynomial. The 
radius of curvature R(z), beam waist w(z) and Gouy 
phase e(z) are given by 



w(z) = two 4/1 + -o-, 



The coefficients C\ and C2 are such that \C%\ <C |Ci|. C2 
depends on the nonlinearity coefficient and length of the 

(2) nonlinear crystal, the magnitude of the pump beam, as 
well as other experimental parameters. The ket |qj, <jj) 
represents a single-photon state in a plane wave mode. 
The vector q, is the transverse component of the wave 
vector kj and <jj is the polarization of the mode j = s 
or i, where the sum is over two orthogonal polariza- 
tion directions Oj and Oj. The polarization state of the 

(3) down-converted photon pair is defined by the coefficients 
C CTs|CTi . The normalized function <I>(q s , q*) is given by *2*3 



R(z) 



and 



e{z) = arctan ■ 



zr 



(4) 
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respectively. The parameter zr is the Rayleigh range. 
The order M of the beam is the sum of the indices N = 
m + n. Note that the usual Gaussian beam is the zeroth- 
order HGoo beam. 

In section IIIII we will make use of the diag onal 
Hermite-Gaussian modes (DHG) defined by ["j****. l20j| 

n+m 

DHG nm (x, y,z) = "5^ b(n, m, k)RGN-k,k(x, y, z), (6) 

with x = (x + y)/ y/2, y — (x — y)/ y2 and the coefficient 
6(n, m, k) defined as 



b(n, m, k) 
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where u(q) is the normalized angular spectrum of the 
pump beam, L is the length of the nonlinear crystal in 
the propagation (z) direction, sinc(a;) = (sinx)/x, and 
K is the magnitude of the pump field wave vector. The 
integration domain D is defined as the region in which the 
paraxial approximation is valid. In most experimental 
conditions, however, D is much larger than the region in 
which $(q s ,qi) is appreciable. 

We assume that 3>(q s , qi) does not depend on the po- 
larizations of the down-converted photons. This assump- 
tion may not be true, especially if the crystal is cut for 
type-II SPDC. However, the polarization dependence can 
be reduced by placing birefringent crystal compensators 
in the down-converted beams |23j . 

We note here that recent experimental work [T21 12*3 . 
*2^ has shown that the quantum state @ is an accurate 
description of the two-photon component of the quantum 
state generated by SPDC using a cw laser. 



II. STATE GENERATED BY SPDC 

Here we review the two-photon quantum state gener- 
ated by SPDC. We consider that a photon p from a suf- 
ficiently weak cw pump beam is incident on a nonlinear 
crystal, producing down-converted signal and idler pho- 
tons s and i, respectively. We will work in the monochro- 
matic approximation, which is justified experimentally by 
the use of narrow bandwidth interference filters in the de- 
tection system. It is assumed that the filters are centered 
at the degenerate wavelength A c = 2A p , where X p is the 
pump beam wavelength. We will also work in the parax- 
ial approximation, which will be discussed below. For a 
sufficiently weak cw laser, the quan tum state generated 
by SPDC can be written as *2lll2!3| 



d I vac) + C 2 



(8) 



where 



'a 3 ,a, JJdqsdqi $(cu,q») \<ls,<J s ) s |q;,0"; 



III. GENERATION OF ENTANGLED HG 
MODES WITH SPDC 

In the following we will denote iv im (q) as the normal- 
ized angular spectrum of the HG mode, which can be 
calculated by taking the two-dimensional Fourier trans- 
form of iJTJ. Explicitly, 

v nm (q x ,q y ) =wD nm H n x 

exp I - A J x 

exp [— i[n + m + l)e(z)] , (11) 
where w = w(z) and 



(12) 



(9) 



such that v nm (q xi q y ) is properly normalized. The general 
problem we are considering is illustrated in FIG.^ We 
now consider that the non-linear crystal is pumped with a 
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FIG. 1: The angular spectrum of the pump beam Vnm(q p ), 
characterized by the wavelength A p and beam width the Wo p , 
creates down-converted fields with angular spectra u,-fc(q s ) 
and v u t{<li), characterized by the wavelength A, 
the beam width wq c = V^wop. 
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Hermite- Gaussian beam HG nm , generating a two-photon 
state \ipnm)- To account for the different wavelengths of 
the pump and down-converted fields, we will write the 
angular spectrum of the HG pump beam as V nm , which 
is equivalent to expression but characterized by the 
wavelength X p and the beam radius Wq p . The angular 
spectrum of the down-converted field v a p is characterized 
by the wavelength A c and the beam radius wq c . It will 
be shown in the appendix that wq c — \Z2wq p . 

Since the HG beams form a complete basis, we can 
expand the two-photon state as 

oo 

\lpnm) = ^ {i(v u t\ s(v 3 k\4>n m )} \Vjk) s \v u t)i , (13) 
j,k,u^t — 

where we have introduced the shorthand notation 



\V a /3, 



= J dqv a!3 (q) |q) 



(14) 



We note here that j (k) and u (t) are the x (y) indices 
of the signal and idler fields, respectively. To facilitate 
the calculations, we will assume that z s — Zi ~ at the 
crystal face. Defining 



we have 



IVv, 



j,k,u,t— 
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(16) 



The task at hand is to calculate the coefficients C^ t . 

For simplicity, we assume that the down-converted 
fields are not entangled in polarization. Then, we can ig- 
nore the polarization dependence of the two-photon state 
(|9"|) . Depending on the type of phase matching, the pump 
and down-converted fields may suffer a small astigmatism 
when propagating through the birefringcnt non-linear 
crystal |25j. This astigmatism depends on the order of 
the modes as well as the length L of the non-linear crys- 
tal, being negligible for thin crystals and/or low-order 
modes. Here we will assume that the crystal is cut for 



type-I phase matching such that the pump beam is polar- 
ized in the extraordinary direction and suffers an astig- 
matism, while the ordinarily polarized down-converted 
fields do not suffer any deformation. We will also assume 
that the pump beam is of low order J\f = n + m < 2 and 
consider that the crystal length is on the order of a few 
millimeters. Under these conditions, we can ignore the 
birefringence and astigmatism effects |25J . Then, using 
Eqs. ©, (HIU and (H3J in Eq. JTBJ gives 



r~i nra 

^jkut 



Vnm (q s sine ( — |q s 



(17) 



In the appendix, we show that 



f t rim 



a\(3\ (I 
At: 



r=0 



arctan A 
(a/2)!(/3/2)! 

-2 

VTTa^ 



b{j,u,a)b(k,t,f3)x 

r 

sinc(r arctan A) 
(18) 



if j + u > n and k + t > m, else CJ%£ t = 0. Here b(j, u, a) 
is given by Eq. Q and we have defined A — L/Kwp], 
a = N — n and [3 — M — m. 

In the appendix, it is also shown that for thin non- 
linear crystals (L <~ 1 mm), Cj%£ t simplifies to 



S * nm 
^jkut 



=\ —6(7, u,N- n)b(k, t, M - m)x 



HGjV-n,M-m(0, 0, 0), 



(19) 
0. Here 



if j + u > n and k + t > m, otherwise Cj%£ t 
HG 7 5 is the Hermite- Gaussian mode Q evaluated at x = 
y = z = 0. When 7 is odd, the Hermite polynomial 
H 7 (0) = 0, which gives another conservation condition: 
N — n and M — m must be even. In other words, the sum 
of the x (y) indices of the down converted fields N = j+u 
(M = k + t) must have the same parity as the x (y) index 
of the pump field n (m). In summary, the conservation 
conditions are 

j + u > n and parity (j + u) = parity n, (20a) 
k + t > m and parity (k + t) = parity m. (20b) 

We note here that the conservation conditions restrict, 
for example, the sum j + u and not j or u individually. 

The parity of the product of the signal and idler HG 
modes (or the sum of the HG mode indices) must main- 
tain the parity of the pump beam angular spectrum, 
which has been transferred to the two-photon quantum 
state. From a mathematical point of view, these results 
are intuitive. For example, consider an even function / 
and an expansion of the sort 



f{x + y) = ^2 Aigi(x)hi(y). 



(21) 
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FIG. 2: Total probability of HG mode generation as a func- 
tion of the order j + k + u + t for a Gaussian pump beams 
with width wo p = 1mm, 0.1mm and 0.05 mm. 



The even parity of / requires that f{x + y) = f(—x — y) 



smaller Wq p , the approximate solution (|19|l is valid only 
for lower orders. Calculations of the total probability 
for an extremely focused pump beam (not shown) shows 
that the total probability for the exact solution converges 
to 1, which indicates that the two-photon state \ip nrn ) is 
properly normalized. 

The parameter of interest is A = L/Kw^ pl which shows 
that the generation efficiency of lower-order modes can 
also be increased by using a longer crystal. However, 
we again emphasize that the pump and down-converted 
fields may suffer greater astigmatic effects in longer crys- 
tals. It is interesting to note that A can also be written 
as A = L/(2zr), where zr is the Rayleigh range of the 
pump beam. This shows that the critical parameter is 
the crystal length L compared to the Rayleigh range of 
the pump beam, arctan^l in (|18|l can be viewed as a 
phase retardation, similar to the Gouy phase J5J. 

Fig. contains the amplitudes C°° Mt and Cjl ut up to 
fourth order (O = 4) for HGoo and HGn pump beams 
with crystal length L = 1 mm and pump beam width 
w p = 0.1mm. For visual clarity, only non-zero terms 
have been included. 



(22) 



Since / is an even function, all products gi(x)hi(y) in 
the expansion must have the same parity, which in this 
example implies that either g and h are both even func- 
tion or g and h are both odd functions. From the point 
of view of physics, the underlying physical process gov- 
erning the generation of HG modes is the transfer of the 
angular spectrum of the pump beam to the two-photon 
state ©, upon which the derivation of the coefficients 
CJkut an d parity and order restrictions above are based. 
We note here that it is also the angular spectrum trans- 
fer which is responsible for the generation of entangled 
orbital angular momentum states 0, 0] . 

We have calculated the exact and approximate proba- 
bility amplitudes for the generation of any arbitrary com- 
bination of HG modes with SPDC. These results show 
that the indices of the HG modes must obey the condi- 
tions PD)|. Eqs. (JT8J, JT5J and are the principal 
results of this paper. Let us now analyze these results 
for some particular HG pump beams with typical exper- 
imental parameters. 

Fig. [21 shows the total probabilities obtained by 



summing all the exact |C*°, ( 

|2 



jkut I 



(circles) or approximate 



\Cjk u t\' 2 (crosses) probabilities up to a given order O 
j + k + u + 1. The pump beam is a Gaussian (n = m = 0) 
with A p = 351 nm and the crystal length L is 1 mm. 
Results are shown for pump beam width wq p = 1 mm, 
0.1mm and 0.05 mm. The total probability £ |C°£ Mt | 2 
approaches unity faster for narrower width pump beams. 
This indicates that experimentally one can increase the 
generation efficiency of lower order modes by focusing the 
pump beam at the plane of the nonlinear crystal. For 



IV. ENTANGLEMENT 

Previous experiments have shown that the pure state 
given by Eq. © is an accurate description of the two- 
photon state generated by SPDC [HHHlil]. Here we 
have show that the two-photon state 10 can also be writ- 
ten as a combination of correlated HG modes in the form 
(|16fl with the normalized coefficients CJ^ t satisfying the 
restrictions on parity and order given in l|20|l . We will 
now use these restrictions to show that the two-photon 
state is entangled in HG modes. 

Let us denote the reduced density operator of - say - 
the signal photon by p^. It is well known that p s has the 
following properties |2a |: (z) p s is a positive operator, 
(ii) tr/3 s = 1 and (Hi) trp 2 < 1 . If p s represents a pure 
state, then trp s = 1, while tip'j. < 1 indicates that p s 
represents a mixed state [2(|. For the pure two-photon 



state 

entangled 2 



trp 2 < 1 implies that the overall state is 
We will show that \ip nm ) is entangled by 



proving that trp < 1. 

It is straightforward to calculate the reduced density 
operator from the two-photon state H16[) : 



Ps = F i kd f \ v jk) ss {v<tf\ , 

jkdf=0 



where 



75=0 



nm 
dfyS- 



(23) 



(24) 



Here we have recognized that the coefficients C™^ t given 
by Eq. I|18|) are rea l- For the argument below, we 
note that: F jkd f = F dfjk , F ]kjk > and trp s = 
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FIG. 3: Coefficients CjlZt up to order O = j + k + u+t = 4 for 
for HGoo and HGn pump beams with width wo P = 0.1mm 
and crystal length L = 1 mm. To improve visualization, only 
non-zero coefficients have been included. 



Srt=o Fjkjk — 1- The reduced density operator p s sat- 
isfies 

oo 

trp 2 s = ( F Jkdf) 2 < 1. (25) 

jkdf=0 

Since p s has unity trace, we can write 

oo oo 

E F 3kjk ^2 Fd f d f = 1 ' ( 26 ) 
jk=0 df=0 



so that from Eqs. 125(1 and (|26|l we obtain 

oo 

E [(^/) 2 -^^M'] <0- (27) 

jfcd/=0 

p s is a positive operator, so its elements satisfy the gener- 
alized Cauchy-Schwartz-Buniakowski inequality |2S| |36j| : 

{Fjkdf) 2 < F jkjk F dfd f. (28) 

Eq. (gU implies that if (F jkdf ) 2 < F jkjk F dfdf for any 
particular values of j, k, d and /, then the equality in 127|) 
must be false, which indicates that trp^ < 1 and \tpnm) is 
entangled. From the parity conservation conditions 121) f> . 
we see that for any 7 in the summation in Eq. I|24|) , j + 7 
and d + 7 must have the same parity as n, otherwise 
Cjk^gCdfjs = 0. A similar relation exists for k + S, f + S 
and to. This implies that Fj kd f — unless j and d and k 
and / have the same parity. The conditions l|2UI) restrict 
the parity of the sum j ' + 7 but not j independently, so 
j can be either even or odd, as is seen in Fig. El for the 
particular cases of HGoo an d HGn pump beams. Then 
there exist Fj k j k ^ and F d f d f ^ such that j and d or 
k and / do not have the same parity, which, using the 
fact that in this case Fj k j k > 0, implies that (Fj kd f ) 2 — 

< FjkjkFdfdf- Then equality in H27JI is false, which 
shows that \ip nm ) is entangled. 

For an infinite dimensional space there will be an infi- 
nite number of terms which satisfy the above conditions. 
The above proof can also be used to show that the state 
\1pnm) is entangled in an arbitrarily large but finite di- 
mensional space, as long as the coefficients (|18|) are prop- 
erly normalized. Experimentally, one can post-select the 
desired HG components of the two-photon state, as will 
be briefly discussed in the next section. As long as the 
reduced density matrix contains one term Fjkdf such that 
j and d or k and / have different parity, the equality in 
(|2Tt is false. 

A. Generating Bell states 

Through post-selection, it is possible to obtain finite- 
dimensional entangled states of higher-order Gaussian 
modes. Experimentally, post-selection can be achieved 
by cou pling th e down-converted fields into optical fibers 

01 HE UH l29| . which filter out unwanted modes. Sim- 
ilarly, entanglement concentration of LG modes was 
achieved by properly coupling these modes into optical 
fibers IH. 

Referring to Fig. for the HGoo pump beam, if one 
considers only first-order down-converted fields (j + k = 
1, u + t = 1), the resulting quantum state is maximally 
ent ang led, resembling the cf> + Bell state, as was observed 
in |16J. It is then fairly straightforward to experimen- 
tally generate all four Bell states using first-order HG 
modes. Using a Dove prism (aligned at 45°) to rotate 
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TABLE I: Amplitudes and probabilities for Hermite- 
Gaussian pump beams HG02 (top) and HG20 (bottom) up 
to second order for crystal length L = 1 mm and beam width 
w p —0.1 mm. The order is defined as O — j + k + u + t. 



HG01 < — > HG10 of either the signal or idler field, one 
can generate the ip + Bell state. Placing one additional 
mirror reflection (or a Dove prism aligned at 90°) in ci- 
ther the signal or idler path, such that HG01 — > — HG01 
and HG10 — ► HG10, one can then generate the maxi- 
mally entangled (\T and ip~ states. 

Another method of generating Bell-states of first-order 
HG modes is with the second-order pump beam HGn. 
Isolating only first-order modes, the output state resem- 
bles the maximally-entangled ip + state, as seen in Fig. 
[3] This method may be advantageous since the high- 
probability zero-order HGoo-HGqo term is not present. 



B. Generating non- maximally entangled states 

Table [I] shows results up to second-order for HG02 
and HG20 beams. Looking at first order terms only 
(HG01 and HG10), the HG02 pump beam creates an 
HG01 — HG01 term, while a HG20 pump generates the 
HG10 — HG10 term. Creating a pump beam that is an 
arbitrary coherent superposition of these two beams, we 
can generate non-maximally entangled pure states. Fig. 
0] shows a possible experimental setup. The input pump 
beam is a horizontally polarized HG02 beam. A half- 
(HWP1) and quarter-wave plate (QWP1) are used to 
adjust the pump polarization. Rotating HWP1, one can 
change the pump polarization as \H) — ► cos#|-ff) + 
sin0|V), where H and V stand for horizontal and ver- 
tical polarization. By tilting QWP1 one can adjust 
the relative phase [30j . The pump polarization is then 
\H) — ► cos9\H) + e^sind\V) up to a global phase. 
The pump beam then enters a polarization-dependent 
Sagnac interferometer with a nested Dove prism (DP) 
orientated at 22.5°. This type of Sagnac interferometer 
is experimentally advantageous since it is insensitive to 
phase fluctuations, and has been used to construct an 
optical single-photon CNOTgate [3l| and to measure 
the spatial Wigner function 32]. It is well known that a 
Dove prism orientated at an azimuthal angle ip rotates an 
image by an angle 2ip in the transverse plane. The polar- 
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PBS ® 45 



HWP1 



1 

□ HWP2 



HG 02 » 

H polarized 



QWP1 



DP @ -22.5 
PBS 



DP @ 22.5 

FIG. 4: Possible experimental setup to generate non- 
maximally entangled states (see text). The input pump beam 
is a vertically polarized HG02 beam. A half- (HWP1) and 
quarter-wave plate (QWP1) are used to adjust the pump po- 
larization. The Dove prisms (DP) rotate the transverse spa- 
tial profile of the pump beam. 



izing beam splitter sends H and ^-polarized components 
into opposite ends of the interferometer, where the Dove 
prism rotates the image of the ii-polarized component 
by 45°, while the ]/-polarized component, which is prop- 
agating in the opposite direction, is rotated by —45°. The 
second Dove prism (located outside the interferometer) 
is used to realign the images in the horizontal- vertical co- 
ordinate system. A Dove prism will also slightly rotate 
the polarization direction. However, since in all cases 
the Dove prisms are followed by PBS's which project 
onto the desired polarization direction, this will result in 
only a slight reduction in beam intensity. After the sec- 
ond Dove prism, the pump beam is in a superposition: 
cos 8 \H) V02 + e 1 ^ sin 9 \ V) V^o- Using a polarizing beam 
splitter (PBS), one can project onto the 45°-polarization 
component, after which the pump beam is in the super- 
position: |45) (cos0Vo2 + e l ^sin0V2o)- The last half- wave 
plate (HWP2) is used to realign the polarization before 
entering the non-linear crystal, where the HG compo- 
nents generate the terms in table [I] Post-selecting the 
first-order terms (j + k = 1 , u + t = 1 ) , the two photon 
state is 



cos 9\v 01 )\v 01 ) + e 10 sin 9 \v 10 ) \v 10 



(29) 



The weights and relative phase of the two-photon state 
(|29|l can be adjusted by rotating HWP1 and tilting 
QWP1, which, along with rotations and reflections of 
the down-converted fields described in the last section, 
allows for the creation of any bipartite pure state. 
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C. Hyperentangled states 



APPENDIX: CALCULATION OF CS 



Another interesting possibility is the creation of hyper- 
entangled states |23| • It has been shown that such states 
may be useful in quantum dense coding and quantum 
cryptography |34j . Using one of the experimental situa- 
tions described above and replacing the single type-I crys- 
tal we have considered with either the type-II "crossed 
cone" source [23J or the two-crystal type-I source [3(| of 
polarization-entangled photons, it should be possible to 
generate a two-photon state entangled in HG-mode and 
polarization. These sources generally require that the 
crystals are thin (on the order of a few millimeters), so 
the possible astigmatism effects discussed in section IIIII 
should be minimal, even for the type-II source. More- 
over, the experimental setups described above require 
only lower order HG modes. 



In section UTTl we showed that the coefficient C^™ t is 



given by Eq. (|T7|> : 



1 2L 

C 3kut =~\hF II dc Ls dc Li v jM<Mx 



V„,„(CK - q, )siiK-( _|q s - qi 



Changing coordinates to 

Q =q s + q i5 

such that dq s dq i = dQ dP/2, we have 



(A.l) 



(A.2) 



V. CONCLUSION 

We have shown that it is possible to generate cor- 
related Hermite-Gaussian modes through spontaneous 
parametric down-conversion. We have derived exact 
and approximate analytical expressions for the probabil- 
ity amplitudes CJ^ t to generate arbitrary combinations 
of Hermite-Gaussian fields. For any Hermite-Gaussian 
pump beam, there exist parity conservation conditions 
for the x and y indices of the down-converted Hermite- 
Gaussian modes. We have used these results to show 
that the two-photon state is indeed entangled in Hermite- 
Gaussian modes. We have discussed the generation of 
maximally entangled Bell-states and non-maximally en- 
tangled pure states of first-order Hermite-Gaussian fields. 
These results can be used to engineer entangled states of 
higher-dimension, and promise to be useful in quantum 
information schemes. 



Q + P\ , /Q P 



V„ m (Q)sinc ( -^P 2 



(A.3) 



Now consider a down-converted HG mode v nm with 
wavelength A c and beam radius wq c . To be more precise, 
let us temporarily write v m „(q, A c , wo c )- Since we are 
working with down-converted fields satisfying A c = 2A p , 
it is easy to show from the general form of HG modes that 
Vnm(q/V2, A c , VZwop) = V„ m (q, X p ,w 0p ). That is, the 
down-converted HG modes with wq c — ^/2wo P will have 
the same Rayleigh range zr, Gouy phase e(z) and radius 
of curvature R(z) as the pump field. Using this property 
of the Gaussian modes, we can expand v* k ((Q + P)/2) 
and v* t ((Q — P)/2) and regroup the x and y terms, which 
gives 



' ik Vut {—) = Vju {-^-' ~v^~) Vkt v^r^ —fr) ' (A - 4) 

where we used definitions @ and l|12|l to show that D* k D^ t = D* u D kt . We note here that relation i|A.4|) is valid for 
all z s — Zi. Then, using the definition of the DHG modes ©, Eq. ljA.41) can be expressed as 



Q f P ' : / Q P ^ = E u ' «) v iV- a ,a (Qx, Px) £ b ( k > 0)VZi-M (Qy> p y) > ( A - 5 ) 

a=0 0=0 



where N — j + u and M = k + t. Noting that D* N _ a a D* M _p a = -D^r-a M-0-^a s> ^ * s straightforward to show that 
the product of angular spectra on the RHS of Eq. (|A.5I) can be rewritten as 

V* N - a , a (Qx,Px) ^M-0,0 (Qy,Py) = V* N _ a>M _ (Q)V* >0 (P). (A.6) 

Putting these Eqs. (|A.5|) and (|A.6|I back into Eq. I|A.3|I . the coefficent Cj%£ t becomes 

i rr N M f f f L \ 

CJ^t = -fcEE^«. Q )HM,|3) / dQV^ M ^(Q)V„ m (Q) / dPV*^(P)sinc — P 2 . (A.7) 

n q=0/3=0 J J \ / 
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The HG modes are orthonormal, so 
J dQ V^ r _ a)M _ | g(Q)V nm (Q) = 5jy- a>n SM-p,m- (A. 8) 
which gives 

C "kut =-\/S?b(3, u,N- n)b(k, t, M - to) 

dP V* N _ nM _ m (V) sine (^P 2 ) , (A.9) 

if N = j + u > n and M = k + t > to, else C$£ t = 0. 

Using the following expression for the Hermite polyno- 
mials 35 1 : 



i/2 



3=0 



(-lpn! 
j!(n-2j)! 



(20 



n-2j 



(A.10) 



it is straightforward to calculate the integral in l|A.9(l 
analytically. After some algebraic manipulation, 



if j + u > n and k + t > to, else CJ^ t = 0. Here we have 
defined A — L/KWp, a — N — n, (3 = M — to and used 
the usual binomial coefficient. 

For thin non-linear crystals, it is possible to arrive at 
a more revealing solution. If the nonlinear crystal is thin 
(L ~ 1 mm), we can approximate smc(L/4KP 2 ) as 1 
in equation (|A~9)l . giving / dPV* i/3 (P)smc(L/4iGP 2 ) « 
JdPV*a(P). Numerical integration shows that errors 
due to this approximation are less than 3% for modes 
as high as a — (3 — 10 for typical experimental values. 
Then, integration in P transforms (|18|l to 



(~l nm 



s~l nm 



-b(j, u,N — n)b(k, t, M — to) x 



HG N—n.M—m 

(0,0,0), 



(A.12) 



if j + u > n and k + t > m, otherwise C, 



jkut 



0. HG 



75 



is the Hcrmite-Gaussian mode Q evaluated at x = y 
z = 0. 



''jkut 



a\/3\ (\ 



Att \2 

(«+/3)/ 2 / a+/3 

E IT 



arctan ^4 
(«/2)!(/3/2) 



■6(j, m, a)b(k, t, P) x 



sinc(r arctan A) 



(A.ll) 
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